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1. Introduction 

If we construct a lattice fermion formulation, there are a number of goals to 
be considered: doubling should be avoided; even at finite lattice spacing a, 
we want to represent chiral symmetry in a sound way; and we are seeking a 
good scaling behavior. Conceptually we have to require locality (the lattice 
Dirac operator D{x^ y, U) has to decay at least exponentially in |x — y| ). In 
addition, for practical purposes we desire a high level of locality, i.e. a fast 
exponential decay or even ultralocality (which means that the couplings 
in D drop to zero beyond a finite number of lattice spacings). A further 
issue is a good approximation to rotation invariance. Last but not least, 
the formulation should be simple enough to allow for efficient simulations. 
Here we report on a construction, which is designed to do justice to all of 
these goals. 

2. Ginsparg- Wilson fermions (an unconventional introduction) 

For a lattice Dirac operator D, full chiral invariance ( {D, 75} = ) is in- 
compatible with other basic requirements (Hermiticity, locality, absence of 
doublers, discrete translation invariance) [1]. Therefore we only implement 
a modified chiral symmetry, which does allow D to fulfill those requirements. 
For such a modified chiral transformation we start from the ansatz 

^ ^ (1 + e[l - Fhs) , (l + e75[l-G])V' , (1) 

e being an infinitesimal transformation parameter. The transformation, and 
therefore F and G should be local, and F,G = 0{a), so that we reproduce 
the full chiral symmetry in the (naive) continuum limit. ^ Invariance of the 



^For convenience, we set a = 1 in the formulae (on an isotropic Euclidean lattice), 
but we classify the terms nevertheless by the order of a that they (would) belong to. 
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Lagrangian il^Dip holds to 0(e) if ^ 

{D,^^} = F^^D + D^^G. (2) 

This imphes a continuous modified chiral symmetry, wliieh has the full 
number of generators. It may be compared to the remnant chiral symmetry 
of staggered fermions: there the doubling problem is not solved, and one 
is only left with a U{\) ® U{\) symmetry, which does, however, protect 
the mass from additive renormalization. The same can be shown here if we 
assume "75-Hermiticity" , = 75-D75, and we choose F = DR, G = RD, 
where R is local again, non-trivial and [i?, 75] = (this generalizes Ref. [2]). 
Then eq. (2) turns into the Ginsparg- Wilson relation (GWR) [3] ^ 

D + D^ = 2D^RD , (3) 

and it implies the absence of additive mass renormalization (see also Ref. 
[5]) since 

{V2RD^V2R-1){V2RDV2R-1) = 1 . (4) 

As another illustration we can write the GWR as {£'"^,75} = 2i?75, and 
we see that a local term R does not shift the poles in (in contrast 
to the cases where {£),75}/2 is local, such as a mass or a Wilson term). 
[With respect to the general ansatz, we have to require the right-hand side 
of D-i + L>t -1 = £)t -ij^Fj^ + GD-^ = D-^F + j5G-f5D^ to be local.] 
75-Hermiticity implies R^ = R. If we now start from some lattice Dirac 
operator Dq (obeying the assumption of the Nielsen-Ninomiya theorem such 
as absence of doublers, but otherwise quite arbitrary), we can construct a 
Ginsparg- Wilson operator D from it by enforcing eq. (4) as 



D 
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A := V2RD0V2R - 1 . (5) 



This is the generalization [7] of the "overlap formula", which uses the "stan- 
dard GW kernel" 

45 := \5.,y , (6) 

^In our short-hand notation, the 'products' are convolutions in c-space. 

^76-Herniiticity is essentially inevitable for any sensible solution, but if we want to 
formulate the GWR even without this assumption, then it reads {D, 75} = 2_D_R75_D. This 
follows from the immediately obvious prescription {Z)~^,75} = 2R'yr,. Alternatively, if 
we require D to be normal we arrive at {Z),75} — 2RD^c^D [4]. However, for the results 
presented in Sec. 3, 4 this doesn't matter, since we always use Rx,y oc dx,y 

*An analogous treatment is also conceivable in the continuum. For comments related 
to the Pauli-Villars regularization, see Ref. [6]. In dimensional regularization 75 is a 
notorious trouble-maker. It may be useful to substitute it by operators (1 — J')76 resp. 
75(1 — G) at the suitable places. In d+e dimensions (d even), F and G could take the form 
eDR/fj, resp. eRD/fj, (where R is some local term, and fj, is the usual scale in dimensional 
regularization). We expect the chiral anomaly to be reproduced correctly as e — » 0. 



3 



and which leads from the Wilson fermion Dq = Dw to the Neuberger 
fermion D = D^e [8]. In the general solution of the GWR, eq. (5), we can 
obviously vary the parameters in Dq or in R (or in both) in many ways, 
without violating their required properties. This shows that there exists a 
continuous set of GWR solutions D in the space of coupling parameters. 

Any solution of the GWR is related to a fully chirally invariant Dirac 
operator = D{1 — RD)~^, which is, however, non-local (in the free case, 
Dy^ip) has poles, cf. eq. (4)). Vice versa, if we start from some with this 
type of non-locality (such as the Rebbi fermion [9], for example) we can 
construct a GW solution [7, 10] D = D^{\ + RD^)~^, which is local, at 
least in the free and weakly interacting case. The mechanism of providing 
locality by inserting a local term i? 7^ is known from the framework of 
perfect actions, where the factor R~^ occurs in a Gaussian block variable 
renormalization group transformation term of the fcrmions [3, 11]. Hence 
— > corresponds to a 5 function block variable transformation, and the 
corresponding perfect action has a Rebbi-type non-locality. The transition 
to locality requires the superficial breaking of the full chiral symmetry, 
R ^ 0: chirality is manifest in the action only in the sense of the GWR, but 
it is fully present in the physical observables [12, 13]. ^ In contrast to the 
Rebbi fermion [16], the axial anomaly is correctly reproduced [17] for the 
perfect action at any local term R, including the perfect Dy. (for R = 0). 
This should also be checked if one generally wants to use D^. in an indirect 
way [18], by measuring the right-hand side of (D^^) = {D^^) — {R). 

By introducing a non-trivial kernel R we have relaxed the condition of 
chiral symmetry somewhat — without doing harm to the physical prop- 
erties related to chirality [3, 5, 2] — and this allows for locality of D (as 
well as the absence of doublers etc.) ^ without contradiction to the Nielsen- 
Ninomiya theorem. In the case of the Neuberger fermion Djsfg, locality 
has been demonstrated in a smooth gauge background. In particular, zero 
eigenvalues in A^A are excluded if the inequality (in d dimensions) 

holds for any plaquette variable P and a suitable bound b. Prom Ref. [19] 
we obtain b = 0.4, which has recently been improved to 6 = {1 + 1/ ^/2)~^ ~ 
0.586 [20]. Still, this constraint is somewhat inconvenient; for instance, at 
least one eigenvalue of A^A has to cross zero if we want to change the 
topological sector. 

^Such a superficial symmetry breaking in tiie transformation term is not necessary in 
order to preserve supersymmetry in a RGT, an lience in a perfect action [14]. A. Tliimm 
also commented on an analogous treatment of further symmetries [15]. 

®This hold at least as long as the gauge background is smooth. At very strong coupling, 
locality is uncertain, and also the doubling problem can return, see Subsec. 4.2. 
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Furthermore, the GWR aUows for locahty only in the sense that the 
couphngs in D{x,y,U) decay exponentially, but not for ultralocality. This 
was conjectured intuitively in Ref. [7]. To demonstrate this No-Go rule for 
GW fermions, it is sufficient to show it for the free fermion. A proof, which 
was specifically restricted to i? = has been given in Ref. [21]. By now, 

a complete proof covering all (local) GW kernels R has been added, hence 
this rule is completely general [22]. 

In that context, it is amusing to reconsider the ordinary Wilson fermion. 
From the mass shift we know that it is certainly not a GW fermion in 
general, and according to our No-Go rule not even the free Wilson fermions 
can obey any GWR. Indeed, if we insert the free Dy[/ into the GWR and 
solve for Rx.y, we find that it decays oc \x — y\~'^ in d = 2, and (x \x — y\~^ 
in d = 4, which is nonlocal and therefore not a GW kernel. 

The exponential decay of D is satisfactory from the conceptual point of 
view, but the existence of couplings over an infinite range is a problem for 
practical purposes. One would hope for a high degree of locality at least, i.e. 
for a fast exponential decay. Almost all the literature on overlap fermions 
solely deals with the Neuberger fermion, but it turns out that the couplings 
in Dj^f, do not decay as fast as one would wish, see Fig. 5. Moreover, 
also other properties listed in Sec. 1 — most importantly scaling, but also 
approximate rotational invariance — are unfortunately rather poor. This 
is obvious even from the free fermion, see Figs. 2, 5. On the level of the 
action, there are generally no 0{a) artifacts for GW fermions, because any 
additional clover term violates the GWR [23], but from the free case we 
see already that the 0{a^) scaling artifacts in Djve are large. ^ And yet its 
simulation is tedious (the quenched case requires already a similar effort as 
simulating Dw with dynamical fermions [25]), allowing only for the use of 
small lattices. 

However, D]\fe arises only from a very special choice in a large class of 
GWR solutions described by eq. (5) , namely i^o = Dw and R = R^^^^ . The 
message of this report is that there are better options, and we are going 
to show in Sec. 3, 4 how im,proved overlap fermions can be constructed, 
tested and applied. This generalization and the improvement concept for 
overlap fermions was introduced in Ref. [7] . It was extensively tested in the 
framework of the Schwinger model [26]. 



'^Our study of the 2d and 4d free fermion, as well as the 2d interaction fermion, show 
consistently that the artifacts in Dns are even worse than those in Dw- Ref. [24] is 
more optimistic about the scaling of the Neuberger fermion (but it does not include a 
comparison to other types of lattice fermions). 
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2.1. THE CONCEPT OF IMPROVING OVERLAP FERMIONS 

We first summarize the main idea: if Dq happens to be a GW operator 
aheady (with respect to a fixed term i?), then eq. (5) yields D = Dq; the 
operator reproduces itself. Therefore, any GW fermion — such as the per- 
fect [3] or the classically perfect [5] fermion — is automatically an overlap 
fermion too. The construction of a classically perfect action for asymptot- 
ically free models only requires minimization — no (numeric) functional 
integral, in contrast to the perfect action — and it still has the additional 
virtues of excellent scaling and rotation invariance. UnfortTinatcly, a really 
powerful quasi-perfect action is not available so far for interacting fermions 
in d = 4. 

However, if we manage to construct at least an approximate GW fermion, 
then we can expect it to change only modestly if we insert it as Dq in the 
overlap formula (for the corresponding R), D Dq. If our approximate 
GW fermion is in addition short-ranged, then we can expect D to have a 
high degree of locality, since the long distance couplings are turned on just a 
little in D. (Also Dq = D]y is short-ranged, but since this is far from a GW 
fermion, it changes a lot in the overlap formula, and those long distance 
couplings cannot be predicted to be tiny). Similarly, if Dq scales well, then 
we can expect this quality to be essentially preserved in D if D « Dq, and 
the same argument applies to the approximate rotation invariance. In Sec. 
3 we discuss examples for promising approximate GW operators. In Sec. 4 
they arc transformed into exact GW fermions, and the above predictions 
are verified. 

3. Short-ranged approximate Ginsparg- Wilson fermions 

Perfect free fermions can be constructed and parameterized in c-space ex- 
plicitly [13]. If we choose the term R such that the locality is optimal at mass 
zero, then we arrive at the standard form R^'^^K To make such a fermion 
tractable, first of all its couplings have to be truncated to a short range. 
A truncation to couplings inside a unit hypercube ( "hypercube fermion", 
HF) was performed in Rcf. [27] by means of periodic boundary conditions. 
(Alternative 4d truncated perfect HFs can be found in Refs. [28].) Trunca- 
tion causes some scaling artifacts, but they are small, so that the free HF is 
still strongly improved over the Wilson fermion. This can be observed from 
the dispersion relation as well as thermodynamic scaling ratios [27, 29]. At 
the same time, truncation also implies a (small) violation of the GWR, in 
agreement with the absence of ultralocal GW fermions. The spectrum of a 
GW fermion with R^^y = 5x,y/{2ii), (/x > 0) is situated on a circle in (Twith 
center and radius fj, (GW circle), see eq. (4). Hence we can test the quality 
of our approximate GW fermion (with respect to i?^**)) by checking how 
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close its spectrum comes to a unit circle. This is shown in Fig. 1 for the 4d 
HF of Ref. [27] on a 20^ lattice. We see that we have a good approximation, 
especially in the physically important regime of eigenvalues close to 0. 




0.5 1 1.5 2 2.5 3 
Re 

Figure 1. The spectrum of a truncated perfect, free HF on a 20* lattice (plotted in W). 

In d = 2 we start off from a similar massless HF, which is optimized 
for its scaling behavior; its set of couplings is given in Ref. [26], and its 
strong improvement over the Wilson fermion is visible in Fig. 2. ^ The free 
spectrum of the "scaling optimal hypercube fermion" (SO-HF) is again close 
to a unit circle, on the same level as the 4d spectrum in Fig. 1. There are 
also other ways to see that we are in the vicinity of a GW fermion, for 
instance by summing over its violation (squared) in each site [26], or by 
inserting the HF into the GWR and solving for R: for the SO-HF, Rx,y 
decays about 6 times faster than the corresponding (pseudo-)i? for Dw- If 
we compare the truncated perfect HF to Dw, then this factor even amounts 
to ~ 71 in d = 2, and to 75 in d = 4. 

We now proceed to the 2-flavor Schwinger model, and we gauge the 
SO-HF by attaching the couplings to the shortest lattice paths only. When 
there are several shortest paths, the coupling is split and attached to them 
in equal parts. Moreover we add a clover term with coefficient 1. For the 

^As an alternative, one could directly optimized the chiral properties by minimizing 
the violation of the GWR [7, 26, 30]. However, those properties can still be corrected 
later on by means of the overlap formula, whereas no such tool is available to correct the 
scaling behavior. Hence it is important to optimize scaling — rather than chirality — 
from the beginning. 

Moreover, if one constructs a truncated approximate GW fermions in a short range, 
there is no direct control over the level of locality of the full GW fermion that one 
approximates. Actually it is not even certain if the GW solution that one would obtain 
by gradually extending the range to infinity is local at all, since also non-local solutions 
of the GWR exist. A simple example for that can be constructed from the SLAG fermion 
[31]. The Dirac operator D = Dslac{^+RDslac)~^ solves the GWR, but it is non-local. 
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Figure 2. Scaling of free fermions: the fermion dispersion relation (left) and the ther- 
modynamic scaling ratio pressure/ [chemical potential)^ at zero temperature (right). 

gauge part we use the standard plaquette action (which is perfect for 2d 
pure U{1) gauge theory [13]). 

This simple "gauging by hand" causes a further deviation from the 
GWR, which is increasingly manifest if the gauge background becomes 
rougher. In Fig. 3 we show the spectra for typical configurations at /3 = 2 
and at /? = 6. It turns out, however, that the SO-HF is indeed an approx- 




0.5 1 1.5 ^ 2 0.5 1 1.5 „ 2 



Figure 3. The spectra of the 2d scaling optimal hypercube fermion (SO-HF) for typical 
configurations at strong resp. weak coupling, approximating a GW circle. 

imate GW fermion up to a considerable couplings strength. Next we test 
the scaling behavior in the presence of gauge interaction. Our simulation 
results (here and below) were obtained in collaboration with I. Hip. They 
are based on 5000 quenched configurations on a 16 x 16 lattice at /9 = 6 
(we use the same set of configurations for all types of fermions), but the 
evaluation does include the fermion determinant, following Ref. [32]. Fig. 
4 shows the dispersion relations for the two meson-type states, a massless 
triplet and one massive mode [33], which we denote as vr and t] (by anal- 
ogy). Again the SO-HF is drastically improved over the Wilson fermion (at 
Kc = 0.25927, from Ref. [34]). It even reaches the same level as a (very 
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mildly truncated) classically perfect action [35], which was parameterized 
by 123 independent couplings per site, whereas only 6 such couplings are 
used for the SO-HF. Therefore, it is realistic to extend the HF formulation 
to QCD, and in fact it has been shown already that minimally gauged 4d 
HFs can indeed be applied in QCD simulations [36, 37]. However, in QCD 
as well as in the Schwinger model, we observed as an unpleasant feature 
of the directly applied HF a strong additive mass renormalization. Using 
the SO-HF, even at /? = 6 the vr mass is renormalized from to 0.13, as 
Fig. 4 shows. This corresponds to a lowest real eigenvalue around 0.03, cf. 
Fig. 3. As an even more striking example, the 4d massless HF, minimally 
gauged by hand and applied to QCD at /5 = 5, leads to a "pion mass" of 
3.0 [27], and at /? = 6 the critical bare HF mass amounts to —0.92 [37] 
(which is unfavorable for the level of locality before truncation, and hence 
for the magnitude of the truncation effects) . We overcome this problem in 
the next section by inserting the SO-HF into the overlap formula. 

[As a further alternative to negative bare mass and overlap, we can reach 
the exact chiral limit solely by the use of fat links [26] . This is essentially 
equivalent to gauging the HF such that the interacting GWR is just violated 
modestly. As yet another possibility one may attach an amplification factor 
> 1 to each link [38].] 




Figure 4- The "meson" dispersion relations for some types of 2d fermions in the 2-flavor 
Schwinger model at P — 6. 

4. Improved overlap fermions 

We now perform the second step in our program and insert the SO-HF 
— which is an approximate GW fermion with an excellent scaling behav- 
ior — into the overlap formula (5) with R = i?^**) given in eq. (6). This 
leads to an exact GW fermion and therefore to all the nice properties re- 
lated to chirality, which are extensively discussed in the recent literature on 
GW fermions: correct anomalies; no renormalization of mass zero, vector 
current and flavor non-singlet axial vector current; no mixing of weak ma- 
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trix elements; no exceptional configurations [5]; correctly reproduced chiral 
symmetry breaking [39] etc. 

Our first prediction was that the level of locality should be improved 
over the Neuberger fermion, and this is clearly confirmed, both, in the free 
and in the interacting case. Fig. 5 compares the decay of the couplings of 
the free fermion (left) and the decay of the "maximal correlation" / over 
a certain distance r — as suggested in Ref. [19] — at /? = 6 (right). (In 
the latter figure one could still try to improve the locality in both cases 
by deviating from and tuning R for optimal locality.) Next we want 




5 10 15 20 5 10 15 20 r 25 



Figure 5. The level of locality for the Neuberger fermions vs. overlap SO-HF: the decay of 
the free couplings in all directions in infinite volume (left) and the decay of the "maximal 

correlation" over a distance r, as defined in Ref, [19], on a 24x24 lattice at = 6 (right). 
The figure on the left also illustrates how well rotation invariance is approximated. 

to verify if the good scaling quality survives the modification due to the 
overlap formula. For the free fermions, this is confirmed in an impressive 
way, see Fig. 2. On fine lattices, the scaling is practically identical before 
and after the use of the overlap formula; only as the lattice becomes really 
coarse, the overlap does some harm to the scaling quality at some point. 

In the presence of gauge interaction, we consider again the "meson" 
dispersions, and we observe again the persistence of the improvement, see 
Fig. 6. 

Furthermore, continuous rotation invariance is approximated much bet- 
ter for the HF than for the Wilson fermion (after all, also this holds exactly 
for the perfect fermion [13]). In Ref. [26] this was also tested in the interact- 
ing case by measuring how smoothly correlations decay with the Euclidean 
distance, and we observed that the SO-HF is again by far superior over 

®One puts a unit source at some site x and defines /(r) := ""y" { ||V'(y)ll | l^: — 2/1 = r }. 
^"The dispersion curves of the overlap fermions just stop inside the Brillouin zone due 
to the square root. The end-points can be shifted by choosing Rx,y = 6x,y/{2jj) and 
varying /u. We show /^ = 1 which is a reasonable choice; if we push the end-point to far, 
the danger of doubling increases. On the other hand, the regime of small momenta is 
needed, of course. 
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Figure 6. The "meson" dispersion relations for two types of overlap fermions at (3 — Q. 

the Wilson fermion; once more it reaches the same level as the classically 
perfect action. The overlap formula (5) suggests that this property is essen- 
tially inherited for the overlap fermions, and indeed we observed a similar 
level for the Wilson fermion and the Neuberger fermion on one hand, and 
for the SO-HF and the overlap SO-HF on the other hand. This confirms 
that also the strongly improved approximate rotational invariance of the 
SO-HF survives if it is turned into an overlap fermion. For the free fermion, 
this progress can also been observed in Fig. 5 (left) from the width of the 



4.1. CHIRAL CORRECTION IN TERMS OF A POWER SERIES 

In d = 2 we can afford an exact evaluation of the notorious square root 
in the overlap formula, but in QCD this is not feasible any more. In the 
recent literature, a number of iterative procedures have been suggested for 
the Neuberger fermion [19, 40]. In Ref. [26] we presented a new method, 
which is very simple and robust, and which is especially designed for the 
case where Dq is an approximate GW fermion already. We evaluate the 
square root in 



the configuration is not extremely rough), hence the expansion converges 
rapidly. On the other hand, for the case of the Neuberger fermion, i.e. for 
Dq = Dw, this expansion fails to converge even in the free case, which 
is presumably the reason why it had not been considered in the earlier 
literature. We call this method a "perturbative chiral correction", where 
the perturbative expansion refers to the GWR violation e (and not to the 
coupling g). It yields a Dirac operator of the form 



'cones 




D„yc = l + AY 



(9) 



11 



For the chiral correction to 0(e"), 1" is a polynomial in A'^A of order n, for 
instance, 

y = [3 - A^A]/2 for n = 1 , (10) 
Y = [15-10A^A + 3{A^Af]/8 for n = 2, etc. 

Hence the computational effort amounts roughly to 1 + 2n matrix-vector 
multiplication (the matrix being A or A^), i.e. it increases only linearly 
(though the convergence is also just linear). Actually this represents a 
fermion with couplings of range 1 + 2n in each component, and it would 
be very tedious to implement it explicitly, even for n = 1. However, due to 
its specific form we never need to do so; -Dp^c can always be evaluated by 
iteration of the above matrix-vector products. 

The crucial question now is if the first few orders are sufficient already 
to do most of the chiral correction. We first look at the free SO-HF, and 
Fig. 7 (left) shows that the first order alone does practically the full job. 
In this context, we also obtain a geometric picture of the effect of the 
overlap formula: for R^^y = 6x,y/{2iJ,) it can be viewed as a projection of 
the eigenvalues onto the circle with center and radius /x. This projection is 
often close to radial. 
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Figure 7. The convergence of the SO-HF to a GW fermion under perturbative chiral 
correction (on a 16 x 16 lattice): the free fermion spectrum (left) and the distribution of 
the small real eigenvalues at p = 6, moving rapidly closer to (right). 

In the interacting case, the convergence to the circle under iteration is 
slowest in the arc around 0. As an example, we show in Fig. 7 (right) a 
histogram of the small real eigenvalues at /? = 6. We see that the mass 
renormalization is removed almost completely if we proceed to 0{£^). 



4.2. BEHAVIOR IN EXTREMELY ROUGH GAUGE CONFIGURATIONS 

For smooth configurations, it is easy to find a parameter — and hence a 
GW circle where the spectrum is mapped on — such that the small (large) 
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real eigenvalues are mapped on This leaves the index unchanged 

(though it is defined by exact zero modes now) and it provides a sensible 
definition of the topological charge via the index theorem [5, 2]. It also 
means that the doubling problem is safely avoided for all typical configu- 
rations at moderate or large (5. 

However, if we dare proceeding to extremely strong coupling, then it is 
not possible any more to find such a center fi oi a GW circle, which does 
the right mapping for all typical configurations. In fact, for extremely rough 
QCD test configurations (on very small lattices) it was observed explicitly 
that all the eigenvalues of the minimally gauged Dhf are close to the arc 
of the GW circle, which is opposite to [41]. For D]^ the eigenvalues are 
densely scattered over a wide area with a large real part. Examples for such 
spectra of an extremely rough configuration are shown in Fig. 8, which was 
provided by N. Eicker, I. Hip and Th. Lippert. At a coupling strength 
where such configurations are frequent, the doubling problem is back for 
those overlap fermions, which are constructed from some simple Dq (for the 
Schwinger model, this problem sets in around /? ~ 1 [41]). This agrees with 
the result of a strong coupling expansion (in the Hamiltonian formulation) 
which applies to D^e [42]. In such cases, the construction of Dn^ maps 
all (almost) real eigenvalues onto the arc close to 2, hence additive mass 
renormalization is back as well — again in agreement with Ref. [42] . In view 
of the latter point one could be tempted to just use a large mass parameter 
/i > 1; this does not help, however, with respect to the doubling problem. 
At very strong coupling, only the (classically) perfect action would help. 

As we mentioned earlier, also locality is in danger in that regime [19], 
and we should therefore keep away from it. As one more advantage of 
choosing Do to be an approximate GW fermion, the regime of /? where we 
are (statistically) on safe grounds is enlarged compared to the Neuberger 
fermion. 

Of course, in the safe regime where (3 is large enough (in QCD this 
includes (3 = 6 for sure [19]) the chiral correction of Dq = Dhf can also 
be carried out by iteration methods different from the one described in the 
previous Subsection, for alternative experiments see Ref. [44] . The efficiency 
in QCD is still to be compared, but for sure in any method the convergence 
will be much faster for Dq = Dhf than for Dq = Dw- 

5. Conclusions 

Our program outlined in Sec. 2.1 has been realized in the Schwinger model, 
and the properties of a resulting improved overlap fermion have been tested 

^^For a Euclidean strong coupling hopping parameter expansion of D^e, see Ref. [43]. 
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Figure 8. The fermion spectrum for an extremely rough QCD test configuration on a 4* 
lattice: for the hypercube fermion of Ref. [27] (left), and for the Wilson fermion (right). 

extensively. They are all clearly superior over the Neuberger fermion, con- 
firming our prediction: the overlap SO-HF scales much better, it is more 
local and it comes much closer to rotation invariance. 

The question now is the applicability of this program in d = 4. The 
4d HF formulation is worked out already, and the corresponding improved 
overlap fermion is currently under investigation in QCD by the SESAM 
collaboration in Jiilich and Wuppertal. 
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